Exact solution of the nonlinear laser passive mode locking transition 
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We present the first statistical mechanics study of a passively mode locked laser which includes 
all the main physical processes, saturable absorption, Kerr nonlinearity, parabolic gain filtering 
and group velocity dispersion, assuming the soliton condition. We achieve an exact solution in the 
thermodynamic limit, where the ratio of the cavity length to the pulse width, the duty cycle, tends 
to infinity. The thermodynamics depends on a single dimensionless parameter 7, the ratio of the 
correlation length to the pulse width. The phase diagram consists of one ordered, mode- locked 
phase and one disordered, continuous wave phase, separated by a first order phase transition at 
7 = 9. The model belongs to a new class of solvable statistical mechanics models with a non-trivial 
phase diagram. The results are obtained with a fully controlled transfer matrix calculation, showing 
rigorously that passive mode locking is a thermodynamic phase transition. 

PACS numbers: 42.55.Ah, 42.65.-k, 05.70.Fh 



Introduction Formation of ultrashort pulses in lasers 
by passive mode locking is an important branch of op- 
tics, both from the point of view of basic research and of 
practical applications. As such, it has been the subject 
of many theoretical and experimental studies over several 
decades, see Q for a review. Nevertheless, until recently 
the central question of the threshold power needed to 
trigger passive mode locking was standing without a sat- 
isfactory answer. 

In several recent works 0, El considerable progress 
was made toward resolving this issue by including the 
effect of noise in a nonperturbative manner. The addi- 
tion of a random clement into the cavity electrodynamics 
turns the laser into a statistical physics system. Apply- 
ing methods of equilibrium statistical mechanics, it has 
been shown that sufficiently strong noise destabilizes the 
pulses formed by a saturable absorber, and that the pro- 
cess of pulse formation is a first order phase transition. 

However, these results were all obtained in the context 
of models of the laser dynamics which assumed a simple 
form of spectral gain filtering. Here, for the first time, 
we tackle the full problem of passively mode-locked laser 
with a fast saturable absorber, quadratic gain filtering, 
slow saturable gain, group velocity dispersion, and Kerr 
nonlinearity, assuming the soliton condition. 

In the first step of the theoretical analysis we construct 
an exact mapping from the statistical steady state of the 
laser electric field to the thermal fluctuations of a string 
with an unstable self-interaction. We show that the ther- 
modynamic limit is obtained when the laser cavity is 
much longer than one of the natural length scales, the 
pulse width or the correlation length. It is next shown 
that thermodynamics is determined solely by the dimen- 
sionless ratio 7 of the two natural length scales. We pro- 
ceed to calculate exactly and explicitly the free energy, 
first by physical, mean field-like arguments, which are 
then established by a rigorous transfer matrix calcula- 
tion. The free energy enables us to calculate all thermo- 



dynamic quantities, including pulse power as a function 
of 7, and the thermodynamic phase diagram, which con- 
sists of one ordered, mode locked phase for 7 > 9, and 
one disordered, non-mode locked phase for 7 < 9. Pulse 
formation is a first order phase transition. The pulsed 
configuration is metastable for 8 < 7 < 9, while the con- 
tinuous wave configuration is metastable for all 7 > 9. 

The model closely resembles equations which have 
been used extensively to study kinetics of phase transi- 
tions and critical phenomena |5(. However, the opposite 
sign of the nonlinearity leads to markedly different phe- 
nomenology. The passive mode locking equation there- 
fore belongs to a new class of statistical physics models. 

Thermodynamics of passively mode locked lasers. Our 
starting point is the master equation which governs the 
slow dynamics of the complex envelope of the electric 
field ip(x,t) in a cavity of length L 
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d t ip = (7 S + ild)d 2 z ip + (7 S -I- i"f k )\^\ 2 ip + gip + rj . (1) 

in < z < L with periodic boundary condition, 
where the real consants j g > 0, 7^, ~f s > 0, 7fc 
are the coefficients of spectral filtering, group velocity 
dispersion, (fast) saturable absorption, and Kerr non- 
linearity, respectively. Noise of spontaneous emission 
and other sources is modelled by the random term 77, 
which is a (complex) Gaussian process with covariance 
(r]*(x,t)r](x',t')) = 2TLS(x - x')5(t - t'). Finally, as 
shown in the slow saturable gain g, may be chosen, 
without significant loss of generality, such that it sets the 
total intracavity power ||V>|| 2 = \ J dx\ip(x)\ 2 to a fixed 
value P. g becomes then a Lagrange multiplier for the 
fixed power constraint. 

In this Letter we consider Eq. in the special but 
important case that 7^ + ijk is a real multiple of 7 9 + ijd, 
known as the soliton condition. In this case one can 
define a 'Hamiltonian' functional IS 



Hty] = £ dx (-~ 7s |</<(z)| 4 + 7 ff l^'(^)| 2 ) , 
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FIG. 1: A typical realization of the real part of the envelope 
i\> of the electric field in the laser cavity, in the mode locked 
phase (top) and in the continuous wave phase (bottom). 



such that the invariant measure p[ip] of equation is 

p[4>] = Z- 1 e- H ^'^5{P~U\\ 2 ) , (3) 
Z= [ [d^][d^*]e- HW/{LT) 8{P- II^H 2 ) . (4) 



Note that the power constraint is enforced explicitly, and 
that p is independent of the imaginary terms in Eq. JQ| . 
The study of steady state preoperties of Eq. is now 
reduced, as in 0, , to the that of an equilibrium statis- 
tical mechanics system with partition function Z. 

The Hamiltonian functional H is almost identical 
to the critical Ginzburg-Landau (GL) functional, the 
paradigm for the effective description of continuous phase 
transition [jj. In Eq. (0), however, in contrast with 
the GL functional, the coefficient of the quartic term is 
negative, making the null configuration unstable. The 
instability is countered by the power constraint, which 
acts nonlocally, and the thermodynamics is consequently 
radically different from the GL one. In particular, the 
one-dimensional system exhibits an ordering transition, 
which is impossible in the one-dimensional GL model. 

The character of the steady state distribution p is de- 
termined by the strength of the ordering saturable ab- 
sorber relative to the disordering noise. The system has 
two natural length scales: The pulse width L p = ^~*p L 
measures the effect of saturable absorption, while the 
correlation length L c — -j^jr measures the effect of noise. 
In both cases, the effect is stronger the smaller is the 
associated length scale. Thus, L p is smaller the larger 
is 7 S and is independent of T, while the converse is 



true for L c . The third length scale in the system, the 
cavity length L, is typically much larger than L p and 
L c in multimode lasers; we therefore study the statis- 
tical problem in the limit L 3> L p , L c , which serves as 
the thermodynamic limit, neglecting small corrections of 
0(L p /L) or 0(L c /L). Thermodynamic quantities, and 
in particular the mode locking threshold are therefore 
determined by the sole dimensionlcss parameter, the ra- 
tio 7 = 4L C /L p = r y s P 2 /T. Note that thermodynamics 
is independent of the spectral filtering 7 S . 

Our analysis proceeds in the textbook approach of cal- 
culating the free energy F = — log Z of the statistical 
mechanics problem Eqs. Hill . from which other thermo- 
dynamic quantities follow. However, seeking to calculate 
the partition function we run into a common obstacle, 
namely that the functional integral in Eq. is not well 
defined. Mathematically this is not a serious problem, 
since the invariant measure is well-defined 0, but in 
order to use Z and F we need to give a precise mean- 
ing to the functional integral, as the continuum limit 
of a regularized version where the integration is finite- 
dimensional. Given a regularization scheme with N ip 
integration, we define the regularized partition function 
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where Hn and || • ||jv are regularized versions of H and 
|| • || . The integration measure is multiplied by factors 
of jjp t° make Zn dimensionlcss, and by factors of ajv, 
a regularization scheme- and iV-dependcnt dimensionlcss 
constant, to make Z — limjv^oo Zn finite. The limit is 
independent of the regularization scheme up to an unim- 
portant multiplicative constant. 

The rest of this Letter is devoted to the calculation 
of F, Eq. and from it the phase diagram and ther- 
modynamic quantities, which all have simple algebraic 
expressions. It is shown that when 7 > 9 the equilib- 
rium is an ordered, mode-locked phase where the power 
P is divided between a single pulse and continuum fluc- 
tuations, see the top panel of Fig. ^ when 7 < 9 the 
equilibrium is a disrodered phase for, where the electric 
field consists only of spatially homogeneous fluctuations, 
see the bottom panel of Fig. Q 

Mean field calculation of the free enegy As a prelim- 
inary step towards the calculation of F we examine the 
problem in two solvable limits. In the first, T — > 0, 
the evolution equation approaches the noiseless dy- 
namics without the random term r\. Then, as is well 
known, ip reaches an equilibrium in the form of a soliton- 

like pulse tp p {z) — e l *y^^sech ^ z ~ z ° ^ , with two initial 

conditions-dependent real parameters, the pulse position 
< zq < L and phase < <j> < 2ir. 

In our context, T — ► is the validity condition for the 
Laplace method, otherwise known as the saddle point ap- 
proximation, in Eq flf. and the solitons ip p take the role 
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H[4, P ]/(LT) 

so that the free energy when L c ^> L p or 7 ^> 1 is 



of the saddle points. It follows that Z ~ e 

T->0 
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In the second solvable limit 7 S — ► 0, becomes 
quadratic and, after replacing the power constraint 
delta function by its Fourier integral representation, the 
functional integral in Eq. (J2J becomes gaussian. The 
quadratic form is then diagonalized using the Fourier rep- 
resentations of ip, and we calculate Z in a regularization 
scheme where ./V is the number of ip Fourier modes kept 
(see Eq. (JjJ). The result of the gaussian integration is 
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Under the assumption that L> L c , we can use the stan- 
dard approximation methods of Euler-MacLaurin sum- 
mation and saddle point integration to get the free energy 
in the limit L c <C L p or 7 -C 1 
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Although the preceding expressions Eqs (l(il-SI) for the 

free energy were obtained in limiting cases, we now argue 
that F p , the pulse free energy, and F c , the continuum free 
energy, may combined into an expression for F valid for 
every 7. The argument is based on the following assump- 
tion: Configurations ip which contribute significantly to 
Z are such that ip(z) = 0(1) for most z, with possibly 
few narrow regions where ip(z) = 0(WL/L p ), whose to- 
tal width is 0(L p ). Let yP, < y < 1 be the total power 
concentrated in regions where ip is large. For z values 
where ip is small the nonlinear term in H is negligible, 
and the existence of regions of large ip affects the statis- 
tics of the small ip region only in that the total available 
power for fluctuations is (1 — y)P rather than P. The 
regions of small ip therefore contribute fLIp^i-yjP to 
the total free energy. Similarly, the regions of large ip are 
so narrow that noise induced fluctuations make negligible 
contribution to the free energy in them, so that the large 
ip regions contribute H[ip]/ (LT) to F. By the principle 
that F is minimized by the Gibbs distribution the large 
ip regions should be such that H[ip] is minimized, that 
is, ip will assume a soliton-like shape with total power 
yP, and contribute F p \p^ y p to F. We claim that since 
the small parameter in these arguments is L p /L, they 
become exact in thermodynamic limit. 

The total free energy is the sum of the pulse and con- 
tinuum parts, minimized over values of y, 
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FIG. 2: The function $(7, y) whose minimum with respect 
to y is the free energy as a function y for 7 = 7, 8, 9, 10 with 
higher values of 7 corresponding to lower curves. Curves with 
two minima correspond to systems with a metastable state. 
When 7 = 9 the values of $(7, ■) at the two minima are equal. 



F can be used to easily derive other thermodynamic 
quantities. An order parameter which is nonzero if and 
only if mode locking occurs is M = y/ (|V , | 4 )- M has di- 
mensions of power, and is proportional to the experimen- 
tally measurable RF power jjj. It follows from Eq. © 
and the definition of F that (IV'I 4 ) = -2T9 7a F. Letting 
$(7,y) denote the target function in Eq. 0, and y(j) 
the minimizer, we calculate 



M( 7 )^-2T(L p /L)^M^ym = (f/3)ip . (10) 



The most important consequence of equation IjlQI) is 
that mode locking occurs whenever the minimizer 1/(7) is 
greater than zero. $(7, •) has the following behavior: For 
7 < 8 it has a single minimum at 0, while for 7 > 8 there 
is a second minimum at i (l+ y/1 — 8/7) , which becomes 
a global minimum when 7 > 9, see Fig- EI Therefore, as 
7 is increased through 7 = 9, M jumps discontinuously 
from to ^-P; that is, the transition is first order, the 
ordered phase is metastable when 8 < 7 < 9, and the 
disordered phase is metastable for all 7 > 9. The pulse 
power in the mode locked phase is y(jy)P, as can be ver- 
ified by calculating higher moment. 

Transfer matrix calculation of the free energy We 
turn to the transfer matrix calculation of Z which estab- 
lishes our expression @ for the free energy with a con- 
trolled derivation. For this purpose we consider a slightly 
generalized partition function Z with fixed values of ip at 
the endpoints of the interval, -0(0) = ipi, ip(L) — ipf, and 
generalized free energy T = — log 2. In this calculation 
it is more convenient to use system parameters that do 
not depend on the system size, so we define 
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Z(iPi,iP f ,a,P,V,L)= [#][#*] x e^o -WOOI'W / da^^a _ 7?) . ( H ) 

Like Z, -Z needs to be defined using a limiting procedure with a properly scaled functional measure. Let Z$ denote 
the lattice regularization of Z with lattice spacing 5. It satisfies the identity 

Zstyi, A V,L)= I -^dWZsWi, 1>, a,[3,V- S\^\ 2 , L - a) e *«M 4 -4l*/-*l a (12) 



We next expand the right-hand-side in powers of S and 
if) — ipf and perform the gaussian integration. This gives 
an equation involving Zg and its derivatives. Taking then 
the continuum limit 6 — > gives the equation for Z 



-d^jd^Z 



\iP f \ 2 &pZ-d L Z + ^a\^ f \ 4 Z = , (13) 



and from it the equation for T 

~|^| 2 - \4, f \ 2 d v T+ ia^/l 4 = , (14) 

where the term with double ipt derivative and the 
term with an L derivative were dropped, since they are 
0(L p /L) with respect to the terms retained. 

Mean field arguments, similar to those presented above 
to obtain F, can be used to calculate T . Like F, T is 
the sum of a continuum contribution T c where the non- 
linear term is unimportant, and a pulse term T v where 
noise is unimportant. Since the continuum fluctuations 
are a bulk property, T c is independent of the boundary 
conditions to leading order, and in this order it is equal 
to F Cy which in the present parametrization reads 



4#P(1 - y) ' 



(15) 



as before, (1 — y) is the relative continuum power. 

Like -Fp, J- p is the negative of the maximal value of the 
exponent in Eq. I|llfl . subject to a given total power yV 
and the boundary conditions. However, the fixed bound- 
ary conditions break the translation invariance, and the 
maximization is achieved when pulses are created near 
the boundaries. We use standard variational methods to 
obtain the result: T v is the minimum of four candidate 
function distinguished by the configuration of the bound- 
ary pulses labelled by the four possible sign choices in 
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^(2Ai±(A-f|^| 2 ) i T(A-f|^| i 



-XyV 
(16) 

In the first ambiguous sign the upper (lower) choice refers 
to the possibility that a pulse maximum lies at z > (z < 
0). The second sign choice refers similarly to the position 
of a pulse maximum relative to the z = L boundary. A 



is a Lagrange multiplier for the power constraint, given 
implicitly by 



yV 



± J\~\a\^\ 2 T\\ 



(17) 

The total free energy T is obtained by minimizing the 
sum T v + T c with respect to A and the possible choices 



of sign, i.e., T = min T c 

A,±,± 



It is left to show that this expression for T satisfies 
the differential equation i|14f> . We demonstrate that each 
of the four possible sign choices satisfies Eq. I|14|) , from 
which the result follows immediately. To this end we 

— j 2 

define A = j- 11 — -t-t, which is the value of A which 

4/37>(l-j/(A)J 

minimizes T for a given choice of sign, as determined 
from Eqs. (|16I17|) The derivatives of T then read 



d v F = \, d^ = Tr f ^l3(X-^f\ 2 ), (18) 

where the upper or lower sign choice corresponds to the 
second ambiguous sign in (|16I17(I . It is now straightfor- 
ward to verify that T and its derivatives satisfy (|14|) for 
all sign choices, as was claimed, finishing the demonstra- 
tion. 
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